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Abstract 

The canard explosion is the change of amplitude and period of a limit cycle born in a 
Hopf bifurcation in a very narrow parameter interval. The phenomenon is well understood 
in singular perturbation problems where a small parameter controls the slow/fast dynam- 
ics. However, canard explosions are also observed in systems where no such parameter is 
present. Here we show how the iterative method of Roussel and Fraser, devised to construct 
regular slow manifolds, can be used to determine a canard point in a general planar system of 
nonlinear ODEs. We demonstrate the method on the van der Pol equation, showing that the 
asymptotics of the method is correct, and on a templator model for a self -replicating system. 
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1 Introduction 



Since the original discovery of canards in the van der Pol equation more than 30 years ago [1], 
they have been identified in numerous systems of nonlinear ODEs. A canard is a trajectory 
which stays close to a repelling slow manifold for an extended amount of time. Canards play 
a key role as parts of transitional limit cycles linking small cycles born in a Hopf bifurcation 
with large relaxation oscillations when a parameter is varied. Since this transition typically takes 
place over a very short parameter interval, and easily may be mistaken for a discontinuous event, 
the phenomenon has been denoted a canard explosion. 

The mathematical theory for canards is well-established for singular perturbation systems of 
the form 

x = f(x,y,c,e), y = eg(x,y,c,e), (1) 

where £ is a small parameter and c is a bifurcation parameter see e.g. [1, 8, 10]. In particular, 
asymptotic expansions in terms of £ of the canard point c c , the parameter values where the longest 
canards exist, can be obtained [13, 3]. However, canard explosions have also been observed in 
many systems that do not have an explicit slow/fast structure with a well-defined small parameter 
e, 

x = F(x,y,c), y = G(x,y,c). (2) 

In some cases a small parameter can be identified after coordinate transformations [ ], while in 
other cases an artificial parameter must be introduced to allow an asymptotic expansion [5, 6]. 
After the expansion, the artificial parameter is set to one to recover the original system. 

While these approaches have been successful, they are somewhat ad-hoc, and it would be 
of interest to establish a systematic approach to identify and locate canard explosions in general 
systems of the form Eqns. (2). The purpose of the present paper is to provide such a procedure. 
It is a simple modification of the iterative method by Fraser and Roussel [ , ] for finding slow 
manifolds. We show that for the van der Pol equation with a distinguished small parameter the 
method gives the correct asymptotic result. For the templator model [ ] with no small parameter 
we get an excellent agreement between the canard point found from simulations and the lowest- 
order canard point from the method. 

2 The canard explosion 

Here we briefly review the basics of the theory for the canard explosion for Eqns. (1). The curve 
defined by f(x,y,c,0) = is denoted the critical manifold S. For £ = 0, S consists of fixed 
points and assuming that it has a fold, the local phase portrait is as shown in Fig. 1. For £ > it 
follows from standard Fenichel theory (see e.g. [ ]) that on the stable side of S an attracting slow 
manifold Ms exists and on the unstable side a repelling slow manifold My exists. The existence 
and the normal hyperbolicity of these manifolds is guaranteed by the theory away from the fold 
point only. However, as trajectories they may be extended across the fold point. In general, Ms 
and Mjj will be distinct, but for a special value of c = c c they may coincide and form a single 
trajectory, a canard. Clearly, the shape of a limit cycle will change dramatically if the parameter 
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Figure 1: Canard explosion in a singular perturbation system, Eqns. (1). As trajectories cross 
the fold of the critical manifold S, they are either repelled down or up, depending on the relative 
position of the slow manifolds Ms and Mu- 



te varied across c c . If M\j is above Ms as in Fig. 1(a) only small limit cycles will be possible. If 
M\j is below Ms as in Fig. 1(b) the limit cycles will be large. 

The single trajectory Ms = My and the corresponding parameter value c c can be found 
asymptotically. For the equation for the trajectories 

f(x,y,c,e)j- = eg(x,y,c,e) (3) 
a Poincare-Lindstedt series is inserted, 

y =yo + yi£+y 2 £ 2 + --- » c c = c + c l £ + c 2 £ 2 + --- . (4) 

Collecting terms of the same order in e algebraic equations for the y^ are obtained. These will 
in general have a singularity at the fold point, but there will be a choice of c^-i such that this 
singularity cancels and yi is well-defined at the fold point. This choice defines the canard point 
and yi is the corresponding canard solution. 

3 A general iterative procedure 

For Eqns. (2) we can also write down the equation for the trajectories, 

F(x,y,c)^ = G(x,y,c). (5) 
Following Fraser and Roussel [9, 1 1], we solve this equation for y algebraically, 

y = *(x,|,c). (6) 

Clearly, it must be assumed that such a solution exists, at least locally. From this equation an 
iterative procedure can be established, 

W = A (7) 



3 



To start the iteration we choose yo such that F(x,yo(x,c),c) = 0, that is, the °°-isocline. Again, 
we must assume that this equation can be solved for yo. Other choices will be possible, but we do 
not have space here to discuss this issue. Typically y^ will have a singularity. Since y# depends 
on c, we will choose the value in each step such that this singularity cancels. This defines the 
procedure for finding canards and canard points for Eqns. (2). 

4 The van der Pol equation 

We now demonstrate the iterative method on the van der Pol equation 

x = y — (x /3 — x), y = e(c—x). (8) 

This system has a canard explosion for c close to 1 when £ is small and positive. The procedure 
from § 2 yields for the canard point [13] 
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cc= l-ie-^e 2 -4^e 3 + ^(e 4 )- (9) 
8 32 1024 



The iterative procedure Eqn. (7) is defined by 

y k+l =x 3 /3-x + e ( —^ (10) 

y k 

with starting point 

y = x 3 /3-x. (11) 

4.1 A numerical example 

We consider first the van der Pol system Eqns. (8) with e = 0. 1 . The asymptotic formula for the 
canard point Eqn. (9) yields c c = 0.986394. 

The iterative process runs as follows: From Eqn. (10) we find 

w^-'+ik^T <12) 

This has a singularity at x = 1 (and also at x = —I, but this is not of interest here) which is 
removed by choosing c = 1, which, then, is the first approximation to the canard point. The 
relative deviation from the asymptotic value is 1.38%. With this choice of c we have 

yi =x 3 /3-x- -L -L- (13) 
\0x+ 1 



and a further iteration yields 



W =*V3-* + (* +1 ^-*> ( 14 ) 

P2 
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where 

p 2 = l(k 4 + 2(k 3 -2(k-9. (15) 

The polynomial p2 has two real roots, x\ = —0.603433 and X2 = 0.987258. We remove the singu- 
larity of y>2 at the latter point by choosing c = X2 = 0.987258, which is the second approximation 
to the canard point. This deviates from the asymptotic value by 0.09%. 
By factorization we get 

p 2 =(x-x 2 )q2 (16) 

where 

q 2 = 10* 3 + 29.8726;t 2 + 29.4919* + 9. 11616 (17) 

such that 

y 2 =x 3 /3-x- { -*^. (18) 

12 

By iteration we find V3, which is a rational function where the denominator is a polynomial of 
degree 8 in x but independent of c. The real roots are -1.24503, -0.999999, -0.389117 and 
X3 = 0.986481. The numerator of V3 is a polynomial of degree 11 in x but linear in c. The 
singularity at X3 can be canceled by choosing c = X3 = 0.98648 1 . This gives yet an improvement 
of the canard point, as the deviation from the asymptotic value is now down to 0.009%. Clearly, 
the procedure can be continued any number of times. 



4.2 Asymptotic analysis 

The structure of the van der Pol equation is sufficiently simple to allow an asymptotic analysis in 
the limit £ — > of the iterative procedure. For a general £, we get in the first iteration 



yi = x 3 /3-x+-^-£. (19) 
x L — 1 



As before, we eliminate the singularity at x = 1 by choosing c = 1 such that 



The next iteration gives 



where 



yi =x 3 /3-x-^—£. (20) 
x+l 



y 2 = x 3 /3-x + {x+1)2{c - x) £ (21) 

P2 



p 2 =x 4 + 2x 3 -2x-l + £. (22) 



The function y2 has a singularity at x = X2 where X2 is a root of P2 and the singularity cancels if 
c = C2 = X2- The polynomial P2 has two roots for £ < 27/16, so the construction of the canard 
trajectory only works under this condition. When £ = 27/16, P2 has a double root at x = 1/2. 
We choose as X2 the root which is greater than 1/2. By a Taylor expansion, one easily finds 

x 2 = c 2 = ^-le--^-e 2 --^-e 3 + 0(e 4 ). (23) 
8 128 2048 v ' 
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This agrees with (9) to 0(e), but not to 0(e 2 ). Inserting c = C2 from Eq. (23) in y 2 and iterating 

in Eq. (10) yields 3^3 as a rational function. In this, we let c = cq + c\£ + C2£ 2 H and in a 

Taylor expansion the first terms are 

cq-x cix 4 + 2cix 3 -(2ci + l)x-c -ci 2 . 

y3 = (x-i)(x+i) £+ (x-mx+ir £+ "' (24) 

By choosing cq = 1 and c\ = —1/8 the singularities at x = 1 in the first two terms cancel. Pro- 
ceeding to the term of order e 3 (we omit the rather long expression) one cancels a singularity by 
choosing c 2 = —3/32. Continuing this way, we find an approximation to the canard point as 

e=l -?-h?-m?+°<rt <25) 

This agrees with (9) to 0(e 2 ), but not to 0(e 3 ). Again, we may continue this procedure to any 
order, in each step correcting a term in the asymptotic expansion of the canard point. 



5 The templator 

The templator is a mathematical model for the kinetics of a self-replicating chemical system. 
The reactions are 

X ^X 
X+X -> T 
X+X + T ^T + T 
T —? P 

The key process the third one where a dimer T acts as a templates and catalyzes its own produc- 
tion from a monomer X. In dimensionless variables the model can be written 



— = r-k u X 2 -k T X 2 T, (26a) 
at 

^ - KX 2 + k T X 2 T - ^— . (26b) 
at K + T 

The last step in the reaction is modeled as an enzymatic reaction with Michaelis-Menten kinetics. 
For further details on the model and its biological significance see [2, 4] and references therein. 

In the following we fix the parameters k u = 0.01, kj = 1, q = 1, K = 0.02 and consider r as 
a bifurcation parameter. In [2, 4] it is shown numerically that two canard explosions occur in the 
model. One is at r = 0.419942 where a small limit cycle explodes as r increases. The large limit 
cycle persists until r = 0.967555, where it turns into a small cycle in another canard explosion. 
See Fig. 2. There is no obvious small parameter in the equations, so the standard asymptotic 
approach for Eqns. (1) does not work. However, in [ ] it is shown that it is possible to account 
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(a) (b) 




Figure 2: Simulations of the templator model Eqns. (26). (a) The black curve is the limit cycle 
for r = 0.419940, the gray curve is is a part of the limit cycle for r = 0.419945. (b) The full large 
limit cycle from panel (a). Note the differences in the scales, (c) The gray curve is the limit cycle 
for r = 0.96755, the black curve the limit cycle for r = 0.96756. 
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for the two canard explosions by two different scalings of the equations. Here we show that the 
iterative method described in this paper can be applied directly on the unsealed equations. 
The equation for the trajectories is 

k u X 2 + k T X 2 T - ^ = r - k u X 2 - k T X 2 T. (27) 

K + T J dT 

This is a quadratic equation in X. Choosing the positive solution, we get for the iteration process 



r(K + T)+qTX' k 
Xk ~ 1 ~ (X^+l)(k u + k T T)(K + TY (28) 

To start the iteration, we choose the dT jdt = isocline as the initial approximation, 



^-^(K + krf^K + TY (29) 

The expression for X\ which is quite complicated has a denominator which is independent of 
r. It has two zeroes, T = 0.0143454 and T = 0.599393. Inserting these in the numerator of 
X\ and requiring that it is zero to cancel the singularity yields a linear equation to determine r 
with solutions r = 0.417681 and r = 0.967710 respectively. The first canard point deviates from 
the numerically determined one by 0.6%, while the latter deviates with 0.02%. Hence, a very 
accurate determination of the canard point is achieved in the very first iteration. 



6 Conclusions 

We have demonstrated that a very simple iteration procedure can be used to determine canard 
points in general planar dynamical systems with no distinguished small parameter. We have 
shown for the van der Pol equation that we obtain an asymptotically correct result in the limit 
of £ —7- 0, and we conjecture this is a general result for problems on the classical singular per- 
turbation form. For the more complex templator system, the method successfully found the two 
canard points in one iteration. In the analytical approach [ ] different scalings were needed to 
find the two canard points. 

It is interesting to note that for the van der Pol equation an upper bound for the small param- 
eter e for canard explosion to occur was found. Recently a bound of £ < 1/4 was found from 
consideration of the curvature of the trajectories [7]. The present bound is more conservative, 
and it would be interesting to obtain a clearer understanding of the relation of the two approaches 
to canard explosion. Furthermore, the present iterative procedure provides a new view on canard 
explosion which may lead to a more general understanding on the specific conditions needed for 
a planar dynamical system without a small parameter to experience a canard explosion. Since 
systems with canard explosions of this kind are abundant in the applications this seems to be of 
fundamental interest. Work along these lines is in progress, and will be reported elsewhere. 



8 



References 



[1] E. Benoit, J. L. Callot, F. Diener, and M. Diener. Chasse au canard. Collectanea Mathe- 
matical 32:37-119, 1981. 

[2] K. M. Beutel and E. Peacock-Lopez. Complex dynamics in a cross-catalytic self-replication 
mechanism. Journal of Chemical Physics, 126:125104, 2007. 

[3] M. Br0ns. Relaxation oscillations and canards in a nonlinear model of discontinuous plastic 
deformation in metals at very low temperatures. Proceedings of the Royal Society of London 
Series A, 461:2289-2302, 2005. 

[4] M. Br0ns. Canard explosion of limit cycles in templator models of self-replication mecha- 
nisms. Journal of Chemical Physics, 134:144105,2011. 

[5] M. Br0ns and K. Bar-Eli. Asymptotic analysis of canards in the EOE equations and the role 
of the inflection line. Proceedings of the Royal Society of London Series A, 445:305-322, 
1994. 

[6] M. Br0ns and J. Sturis. Explosion of limit cycles and chaotic waves in a simple nonlinear 
chemical system. Physical Review E, 64:026209, 2001. 

[7] M. Desroches and M. R: Jeffrey. Canards and curvature: the 'smallness' of e in slow-fast 
dynamics. Proceedings of the Royal Society of London Series A, 467:2404-2421, 201 1. 

[8] W. Eckhaus. Relaxation oscillations includig a standard chase on french ducks. In Asymp- 
totic Analysis II, volume 985 of Lecture Notes in Mathematics, pages 449-494. Springer 
Verlag, New York/Berlin, 1983. 

[9] S. J. Fraser. The steady state and equilibrium approximations: A geometrical picture. Jour- 
nal of Chemical Physics, 88(8):4732-4738, 1988. 

[10] M. Krupa and P. Szmolyan. Relaxation oscillations and canard explosion. Journal of 
Differential Equations, 174:312-368, 2001. 

[11] M. R. Roussel and S. J. Fraser. Geometry of the steady-state approximation: Perturbation 
and accelerated convergence methods. Journal of Chemical Physics, 93(2): 1072-1081, 
1990. 

[12] F. Verhulst. Methods and Applications of Singular Perturbations. Number 50 in Texts in 
Applied Mathematics. Springer, New York, 2005. 

[13] A. K. Zvonkin and M. A. Shubin. Non-standard analysis and singular perturbations of 
ordinary differential equations. Russian Mathematical Surveys, 39(2):77-127, 1984. 



9 



